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Abstract 

The lifetime difference (ud = AT £>/2T d) and mass difference {xe> = Amn /T jj) of neutral D meson 
have been measured with yo = (0.80 ±0.13)% and xd = (0.59 ± 0.20)%, respectively. Intriguingly, 
in contrast with the cases of K and B q systems, the current data indicate that vd/xd ~ 1 and 
yu favors to be larger than xq. For explaining the experimental indication, we here study the 
D — D oscillation in the framework of unparticle physics. We demonstrate that the peculiar phase 
appearing in off-shell unparticle propagator could play an important role on xd and yp. 
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In the Standard Model (SM), the most impressive features of flavor physics are the 



Glashow- Iliopoulos-Maiani (GIM) mechanism [l| and the large top quark mass. The former 
results in the cancelation between the first two generations so that the mass difference Am^ 
in the neutral K system could be suppressed, while the latter makes Arris (q = d, s) in the 
B q systems dominated by the short-distance (SD) top-quark effects j^j. Due to the precision 
measurements and the sensitivity to the new physics, within the past decades enormous 
studies have been done in K and B q mesons, which are composed of down type quarks. 

By the production of large number of D mesons at Tevatron and B factories worldwide, 
now the neutral charmed meson which is made up of up-type quarks also plays an important 
role on the test of the SM. By the world average, the current measurements with allowing 
CP violation (CPV) for D — D mixing are given by jsj 

J- D J- D 

Vd = = =£ = (0-80 ± 0.13)% . (1) 

Combining the errors in quadrature, the ratio of yu to xd is estimated by 

— = 1.36 ±0.42. (2) 

x D 

Intriguingly, the current data not only show ~ yr> but also indicate that the former is 
slightly smaller than the latter. 

Due to the effective GIM mechanism and the absence of heavy quark enhancement, the 
SD SM predictions are several orders smaller than the data j^]. It is expected that the GIM 
suppression factor might be lifted by long-distance (LD) effects 0-8]. With the exclusive 
technique EJ, the results in the SM are estimated to be js] 

x D M « (0.108 ± 0.05)% , y D M « (0.30 ± 0.34)% , (3) 

where we have averaged the possible theoretical scenarios. Since the exclusive technique 
is based on the measurements of nonleptonic D decays, due to the limited accuracy of 
experimental data, the SM prediction on yu is still quite uncertain. Although the results in 
Eq. ([3]) display the same tendency as the data, the values of and y|, M are quite smaller 
than the experimental data. Thus, the ratio in the SM is estimated as 

ySM 

- 2.78 ± 3.40 . (4) 
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We see clearly that the central value by LD contributions is twice larger than that in Eq. ((2)). 

If we take the central values of data in Eq. ([I]) seriously, the SM results in Eqs. fl3]) and 
fll]) obviously cannot match with the data consistently. For explaining the large xo{vd) 
and Dd/%d ~ 1, the incompatibility could be ascribed to new physics. In most extensions 
of the SM, owing to the suppression of (m c /my/) 2 j^] and the constraints of low-energy 



measurements 
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11] . the SD contributions to yo with O(10 3 ) is not favorable. Therefore, 



we are going to explore a peculiar new effect on the D — D mixing, especially on the yp, 
where the associated new stuff is dictated by the scale or conformal invariance and named as 



unparticle 
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referred to Refs. 



13 L Some interesting applications of unparticle to various systems could be 
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18] . The unique character of unparticle is its 



peculiar phase appearing in 



the off-shell propagator with positive squared transfer momentum 



12]. Due to CP invariance, 



the ima gina ry (real) part of the phase factor leads to the absorptive (dispersive) effect of a 



process 



19 



20] . In this Letter, we investigate how xr, and yp are influenced by the phase 



factor. Furthermore, in order to make the production of scale invariant stuff be efficient at 
Large Hadron Collider (LHC ), w e will concentrate on the unparticle that carries the color 
charges of SU(3) C symmetry [17]. 

Since there is no well established approach to give a full theory for unparticle interactions, 
we study the topic from the phenomenological viewpoint. In order to avoid fine-tuning the 
parameters for flavor changing neutral currents (FCNCs) at tree level, we assume that 
the unparticle only couples to the third generation of quarks before electroweak symmetry 
breaking. Hence, the interactions obeying the SM gauge symmetry are expressed by 

^7 teWOJ + IlQlI^Q^O^] , (5) 

where Ir^l are dimensionless free parameters, q' R = tR, 6r, Q t l = (t,b)i, {T a } = {X a /2} are 
the SU(3) C generators (where A a are the Gell-Mann matrices) normalized by tr(T a T b ) = 
S ab /2. Kjj is the scale below which the unparticle is formed, and the power djj is determined 
from the effective interaction of Eq. (jSJ) in four-dimensional space-time when the dimension 
of the colored unparticle Oy is taken as dy. Since we only concentrate on the phenomena 
of up type quarks, the associated interactions are formulated by 

U lfi (XrPr + X L P L ) T a Ud^O v , (6) 

where U T = (u,c,t), ^r(l) is a 3 x 3 diagonal matrix and diag(X fi( x)) = (0, 0, Ir^/A^). 
After spontaneous symmetry breaking of electroweak symmetry, we need to introduce two 



unitary matrices V v ' to diagonalize the mass matrix of up type quarks. In terms of physical 
eigenstates and using the equations of motion, the interactions for c—u — Oy could be written 

as 



C 



m 



c - I R 



(9ucPl + 9 L uc Pr) T a cO a v + h.c. 



(7) 



where the mass of light quark has been neglected. And g* c = A x ( V(j )i3(V^ ) 2 3 with x = R,L, 
in which the index of Arabic numeral (1, 2, 3) stands for (u,c,t) quark, respectively. 



By following the scheme shown in Ref. 
is written as 



18j | , the propagator of the colored scalar unparticle 



d* xe - lkx (0\TO a (x)O b (0)\0) = i- 



C s S ab 



-k 2 - ie) 2 ~ du 



(8) 



with 



A 



da 



2 sin dyn 

16tt 5 / 2 T(du + 1/2) 



(9) 



(2vr) 2 ^ F(du - l)T{2du) ' 
Combining Eqs. ((7|) and ([8]), the four fermion interaction for D-meson oscillation is given by 



C s (mi 



2ml V Al 



2 \ d u 



-id\j~K 



[u {g« c P L + gLPR)T a c}- 



(10) 



For estimating the transition matrix elements, we use 

(D\uP R{L) cuP R{L) c\D) 
(D\uP R cuP L c\D) 

(D\u a P R C,3U(3P L C a \D) 

(D\u a P R{L) c p upP R{L) c a \D) 
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1 1 

24 + 4 
1 1 
8 + \2 



^Dm D fl , 



777 + -tZd ) m D fl , 



^ + ) m D f D , 



24 



t,Dm D f D 



(11) 



where £d = m 2 D /{m c + m u ) 2 and fo is the decay constant of D meson. As a consequence, 
the dispersive and absorptive parts of D — D oscillation in the unparticle physics are found 
by 



#12 = M 12 ~ 2 r< L2 > 
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where M^ 2 = cos(duir)hu and Ij 2 = 2sm(du'^)hu with 

36m 2 c \Af 

In order to study the xd and yp, we have to know their relations to M 12 and T 12 . Following 



C s [ml 



(12) 



the notation in Particle Data Group (PDG) 2l|, the mass and rate differences of heavy and 
light D mesons could be formulated by 

Am D = Re(Au HL ) , 
AT D = -2Im{Au HL ) (13) 

with 




Aco HL = 2 X I[M 12 - t -T 12 j {M* 12 - l -T* 12 j , (14) 

where M 12 = Mf 2 M + and T 12 = Tff + Tf 2 . If we define the relative phase between M 12 
and Ti 2 to be <f>D = arg(Mi 2 /ri 2 ), the ratio of rate difference to mass difference is obtained 
by 

AT D 2r D 



Am D 1 - r 2 D /A + R D 



COS (15) 



with 



12 



M- 



12| 



Rd = ^(l-rl/4) 2 +r 2 D cos<f> D . (16) 

We note that unlike the case in B q system where the sign of AT s q in the SM is certain 
and experimental data are consistent with SM prediction, the sign of AI^ in the SM is 
uncertain; thus we use <pD = arg^M^/Tw) for D-meson, instead of 4>b = cirg(— Mi 2 /rf 2 ) 
for Sg-meson. Hence, the ratio of yo to xd can be expressed by 

Up _ Afj, r D cos <p D 

x D ~ 2Am D ~ 1 -r 2 D /A + R D ' 1 ' 

In order to illustrate the phase effect of unparticle and simplify the numerical estimates, 
we set Ajy = 1 TeV and g^ c = g^ c = \9uc\z %b \ the couplings are vector-like. Since the SM 
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predictions are still quite uncertain, for numerical analysis we adopt the recent SM results 
to be |s | 



M- 



SM 
12 

12 



0.13% ps^ 1 
0.73% ps~ x 



(18) 



where we adopt Mf^ 1 = x^To/2 and Ff^ 1 = y^Tm and we take only the central value of 
x i) M (l/f) M ) as input. Other relevant values used for numerical estimates are listed in Table fl] 



TABLE I: Values used for numerical estimates [21]. 



m D [GeV] 


m c [GeV] 


Sd [MeV] 


T D [PS] 


1.864 


1.3 


206.7 


0.41 



With the chosen scenario for the free parameters and the taken numerical values, now 
we have to deal with three free parameters, i.e. the scale dimension du, the magnitude of 
coupling g uc and its phase 9. Since the SM results are smaller than the current data, we 
find that the influence of 9 is insignificant when the constraints of measured xd and yo are 
included. In Fig. [TJ we present the unparticle contributions to and yo as a function \g uc \ 
(in units of 10~ 2 ) and du, where figure (a)-(d) stands for 9 = (0, 7r/4, tc/2, 37r/4) and solid 
and dotted line denotes xd and yu, respectively. It is clear that the allowed \g uc \ is slightly 
changed when 9 is varied. For further understanding the ^-dependence, we plot x& and yu> 
as a function of 9 and du with \g uc \ = 1.5 x 10~ 2 in Fig. [2} where the solid and dotted line 
corresponds to and yo, respectively. 

We have studied the mixing parameter and lifetime difference of D — D oscillation in 
the framework of unparticle physics, where the new stuff is dictated by scale or conformal 
invariance. Unlike other models, due to the peculiar phase of unparticle, not only the mix- 
ing parameter X£> but also the lifetime difference yo can be enhanced to fit the current 
experimental data, especially the experimental result of yo/xD ~ 1- We speculate that the 
unparticle or unparticle-like effects could be strongly verified, when xd ~ yn ~ few x 10 -3 
and yo/xD ~ 1 are satisfied simultaneously in experiments. 
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FIG. 1: (a)-(d) the contours for Amp (blue solid) and AI^ (red dotted) as a function \g uc \ and 
djj with 9 = 0, 7r/4, 7r/2, 4ir/4, respectively. The numbers on the curves are the data with la errors. 




1.45 



FIG. 2: xd (blue solid) and yjj (red dotted) as a function of 6 and djj with \g uc \ = 1.5 x 10~ 2 . 
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